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Анотація. Задачі оптимізації розміщення об'єктів мають широкий спектр застосування. Одним з таких 
застосувань є задача розміщення контейнерів з відпрацьованим ядерним паливом (ВЯП) на площадці зберіган- 
ня. Розв'язання задачі може бути зведене до розв'язання задачі пошуку оптимального розміщення заданого на- 
бору конгруентних кіл у багатозв'язній області з урахуванням заданих технологічних обмежень. У роботі по- 
будовано математичну модель задачі та розглянуто її особливості. Запропонований підхід заснований на мате- 
матичному моделюванні відносин між геометричними об'єктами за допомогою методу рбі-функцій. Це дозво- 
лило звести вирішення задачі до задачі нелінійного програмування. На сьогоднішній день важливою науковою 
проблемою є проблема створення умов для безпечного зберігання відпрацьованого ядерного палива. У про- 
цесі створення будь-якого сухого сховища відпрацьованого ядерного палива можна виділити наступні основні 
етапи: вибір майданчика розташування, проєктування сховища, будівництво, експлуатація та зняття з експлуа- 
тації. Повна перевірка на відповідність сховища та його елементів зазначеним нормам зазвичай розпочинає 
здійснюватися на етапі проєктування. На етапі ж вибору майданчика розташування перевірка на дотримання 
норм безпеки здійснюється лише з точки зору впливу сховища у цілому на навколишнє середовище. Такий 
підхід не можна вважати у повній мірі доцільним, адже, врахувавши, наприклад, усі кліматичні особливості те- 
риторії розташування майбутнього сховища, можна корегувати теплові режими зберігання відпрацьованого 
ядерного палива. Аналогічно можна вважати необхідним аналіз та вибір форми майданчика зберігання з метою 
розміщення максимально можливої кількості контейнерів з ВЯП. Такий вибір, вочевидь, повинен здійснюва- 
тися із врахуванням норм ядерної, радіаційної та теплової безпеки, а також із дотриманням технологічних об- 
межень. Проблему пошуку оптимального розміщення контейнерів з урахуванням заданих технологічних обме- 
жень можна сформулювати у вигляді задачі оптимізації геометричного проєктування. Тому метою досліджен- 
ня є побудова математичної моделі задачі та дослідження її характеристик для розробки ефективних методів 
розв'язання. Запропонований підхід базується на математичному моделюванні відносин між геометричними 
об'єктами за допомогою методу рбі-функцій. Це дозволило звести вирішення задачі до задачі нелінійного 
програмування. 

Ключові слова: відпрацьоване ядерне паливо; ядерна та теплова безпека; задачі оптимізації упаковки; 
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Вступ 

На сьогоднішній день важливою нау- 
ковою проблемою є проблема створення 
умов для безпечного зберігання відпрацьо- 
ваного ядерного палива (ВЯП). Найпоши- 
ренішою та економічно обгрунтованою є 
технологія так званого «сухого» зберіган- 
ня, яка передбачає розміщення відпрацьо- 
ваних тепловиділяючих збірок у спеціаль- 
них ємностях, наповнених інертним газом, 
та подальше їх зберігання у відкритих або 
закритих зонах зберігання. Створене схо- 
вище повинно відповідати всім стандартам 
безпеки, прийнятим для об'єктів цього ти- 
пу. Зокрема, не перевищення загального 
рівня іонізуючого випромінювання та, як 
показують останні дослідження, дотриман- 
ня режимів теплового зберігання відпра- 
цьованого ядерного палива. 

У процесі створення будь-якого сухо- 
го сховища відпрацьованого ядерного па- 
лива можна виділити наступні основні ета- 
пи: вибір місця зберігання, дизайн схови- 
ща, будівництво, експлуатація та виведен- 
ня з експлуатації. Повна перевірка відпо- 
відності сховища та його елементів цим 
стандартам зазвичай починається на стадії 
проєктування. На етапі вибору місця збері- 
гання перевірка на відповідність стандар- 
там безпеки проводиться лише з точки 30- 
ру впливу сховища в цілому на навколиш- 
нє середовище. Цей підхід не можна вва- 
жати повністю доцільним, оскільки, врахо- 
вуючи, наприклад, усі кліматичні особли- 
вості розташування майбутнього сховища, 
можна регулювати режими теплового збе- 
рігання відпрацьованого ядерного палива. 
Подібним чином, можна вважати необхід- 
ним проаналізувати та вибрати форму пло- 
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щадки зберігання, щоб вмістити макси- 
мально можливу кількість контейнерів з 
ВЯП. Такий вибір, очевидно, слід робити з 
урахуванням норм ядерної, радіаційної та 
теплової безпеки, а також з дотриманням 
технологічних обмежень. 

Проблему пошуку оптимального роз- 
міщення контейнерів з урахуванням зада- 
них технологічних обмежень можна сфор- 
мулювати у вигляді задачі оптимізації гео- 
метричного проєктування. 

Теорія оптимізаційного геометрично- 
го проєктування є однією із інструментів 
вивчення та оптимізації складних техніч- 
них систем з метою досягнення стану їх 
оптимального функціонування. Вона при- 
значена для розв'язання ряду прикладних 
оптимізаційних задач розміщення геомет- 
ричних об'єктів. Ці задачі пов'язані зі ство- 
ренням енерго- та ресурсозберігаючих тех- 
нологій в пріоритетних галузях народного 
господарства (енергетиці, машино-, судно-, 
авіабудуванні, будівництві, хімічній про- 
мисловості, а також у наукових досліджен- 
нях в області нанотехнологій, в сучасних 
задачах біології, мінералогії, медицини, ма- 
теріалознавства, у робототехніці, задачах 
кодування інформації, системах розпізна- 
вання образів, системах керування косміч- 
ними апаратами) при автоматизації та 
моделюванні процесів розміщення різних 
об'єктів. 

Метою дослідження є побудова 
адекватної математичної моделі задачі по- 
шуку оптимального розміщення контейне- 
рів з ВЯП та дослідження її основних влас- 
тивостей, які визначать напрями розробки 
ефективних методів пошуку розв'язку. 
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Аналіз останніх досліджень |і 
публікацій 

Проблема зберігання відпрацьовано- 
го ядерного палива (ВЯП) у світі, зокрема в 
Україні, останнім часом стає все більш ак- 
туальною (1-2). Оскільки частка ядерної 
енергії в загальній кількості виробленої 
електроенергії традиційно висока й переви- 
щує 50 95, щороку енергетичні реактори 
чотирьох діючих атомних електростанцій 
(АЕС) утворюють значну кількість високо- 
активних радіоактивних відходів та ВЯП. 
Проблема поводження з ВЯП енергетич- 
них реакторів в Україні регулюється, 30- 
крема, наказом Міністерства енергетики 
України від 03 вересня 2020 р. Мо 566 «Про 
затвердження Стратегічних напрямів по- 
водження з відпрацьованим ядерним пали- 
вом атомних електростанцій України з ре- 
акторами типу ВВЕР на період до 2030 
року та Плану заходів щодо їх реалізації». 
Він передбачає тривале зберігання ВЯП на 
ВП Запорізька АЕС (ЗАЕС), паливо інших 
діючих атомних електростанцій -- у центра- 
лізованому сховищі ВЯП Чорнобильської 
зони відчуження. Тривале зберігання ВЯП 
на території України сьогодні фактично 
здійснюється лише на найбільшій АЕС - 
Запорізькій. Пристанційне сховище ЗАЕС 
призначено для зберігання понад 9 тисяч 
відпрацьованих  тепловиділяючих збірок 
шести реакторів ВВЕР-1000 та розраховане 
на експлуатацію близько 50 років, тобто 
протягом передбачуваного терміну експлу- 
атації ЗАЕС. Окрім ВЯП ЗАЕС, в Україні 
планується ввести в експлуатацію ще два 
сховища ВЯП. Це сховище для палива ре- 
акторів РБМК-1000 Чорнобильської ДЕС 
(СВЯП-2) та централізоване сховище 
(ЦСВЯП) для палива інших АЕС України. 
Обидва ці сухі сховища застосовуватимуть 
так званий "сухий" спосіб зберігання, тобто 
без використання води в якості охолоджу- 
вача. Метод сухого зберігання ВЯП є до- 
сить поширеною стратегією поводження з 
відпрацьованим паливом для країн з від- 
критим ядерним паливним циклом. Існує 
декілька типів сухих сховищ, два з яких бу- 
дуть використовуватися в Україні - кон- 
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тейнерне сховище на ЦСВЯП та модульне 
(на СВЯП-2). Незважаючи на різницю в 
конструкції основного обладнання для 
зберігання, сухі сховища мають спільні 
науково-технічні питання безпеки експлуа- 
тації, на які слід звертати увагу та спрямо- 
вувати зусилля на їх наукову підтримку. 

Серед таких проблем можна виділити 
проблему оптимального використання пло- 
щі зберігання. Вирішення проблеми пошу- 
ку оптимального розміщення контейнерів з 
ВЯП може бути зведене до вирішення 
проблеми пошуку оптимального розміщен- 
ня заданого набору конгруентних кіл у ба- 
гатозв'язній області з урахуванням техно- 
логічних обмежень. Сьогодні багато дос- 
ліджень присвячені задачам упаковки рів- 
них кіл у квадрат. Проблема вивчалась у 
великій кількості статей як з теоретичної, 
так 1 з обчислювальної точки зору. До- 
слідження включають алгоритми промене- 
вого пошуку |3|, автоматизовані докази оп- 
тимальності |4| та підходи, що використо- 
вують методи нелінійного програмування 
51. 

Різні евристичні підходи для роз- 
в'язання задачі розміщення рівних кіл у 
круглому контейнері були запропоновані 
в (6, 71. 

В даний час багато статей розгляда- 
ють задачу упаковки фіксованого набору 
об'єктів у заданий об'єкт, з метою мініміза- 
ції його розмірів. Кілька статей вирішують 
цю проблему за допомогою нелінійних мо- 
делей та методів оптимізації. У роботі |З| 
наведено огляд таких досліджень. Також у 
|8) розглядається велика різноманітність 
моделей нелінійного програмування для 
задачі мінімізації розмірів об'єктів у 2Д- 
задачах упаковки (квадрат, прямокутник, 
трикутник та коло). 

Слід відзначити, що відомі нам до- 
слідження не вирішують задачу упаковки 
рівних кругів у багатозв'язну область, ме- 
жа якої складається з дуг кіл 1 відрізків пря- 
мих ліній. 
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Постановка проблеми 

Загальну постановку оптимізаційної 
задачі геометричного проєктування можна 
сформулювати таким чином: визначити, З 
урахуванням всіх правил, вимог та техно- 
логічних обмежень таке просторове розмі- 
щення визначеного набору геометричних 
об'єктів в заданому контейнері, при якому 
критерій оптимізації приймав би екстре- 
мальне значення. 

Таким чином, для постановки задачі 
розміщення контейнерів з ВЯП в термінах 
задачі оптимізаційного геометричного про- 
єктування необхідно в аналітичному вигля- 
ді визначити: 

1. просторову форму області розміщення; 

2. просторову форму об'єктів, які необ- 
хідно розмістити в області розміщення; 

3. технологічні умови розміщення зада- 
ного набору об'єтів у межах зони роз- 
міщення; 

4. критерій оптимізації розміщення кон- 
тейнерів. 

Формалізуємо кожен з перелічених 
пунктів. 

Просторова форма області роз- 
міщення. Зовнішній контур геометричної 
форми Р може бути визначений у відпо- 
відності до топографічних умов обраної 
місцевості. 


Рис.1. Врахування топографічних умов 
обраної місцевості при проєктуванні 
майданчика 
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Топологічна форма області Р може 
бути однозв'язна або багатозв'язна. Май- 
данчик Р може бути багатозв'язною об- 
ластю через те, що в середині майданчика 
можуть знаходитись участки, в яких не- 
можливо розмістити контейнери (наприк- 
лад, це може бути службова зона обслуго- 
вування чи особливості рельєфу). 

Для математичного моделювання за- 
дачі аппроксимуємо зовнішню межу май- 
данчика набором відрізків прямих ліній та 
дугами кіл. Припустимо, що Р формуєть- 
ся наступним чином: 


С 
Ресі| ВМ )А, |; 
І-І 
де сі() є замикання (); В; являє собою 


просту зв'язану множину, межа якої 
утворена послідовністю відрізків лінії 


Ї ху» УреХни» Ура | і дуг кіл 
| вУваавУньхрУрТі Р де (х,У;) і 
(уз Урі) Є координатами початку та 
кінця відрізка або дуги, (Х;, У;) і Т; Є 
центральними координатами та радіусом 
кола відповідно; А; є зоною заборони на 


розміщення контейнерів та може бути 
представлена як 


1 У, 
ла а, Ш Ом 
2-1 4-1 


рей а1Б201) 
Сі - (ху) є Е? ох) КОЗУ) з 


(ру) «05,8 єї, є (2, 


Мід - опуклий багатокутник, заданий 


Ч 
тд вершинами, 4 є Із - П,2,...,3), тобто 
кожна зона заборони - це не опукла мно- 
жина, яка може бути представлена кінце- 
вим об'єднанням різних кіл та опуклих 
багатокутників (рис. 2). 
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(ХУ) ЛЯ 
ль ьЗ (ху ) (ХУ; ) 


(х,У,), 


(ХУ, 


(Хо) 


б 


Рис. 2. Приклад вигляду області Р 
а) однозв'язна область; 
6) багатозв'язна область 


Описаний спосіб задання області Р 
дозволить з певним рівнем наближення 
описати будь-яку форму. 

Просторова форма об'єктів, які не- 
обхідно розмістити в області розміщення. 
Нехай задано набір з М циліндричних кон- 
тейнерів К., і є Гу , З ВЯП та майданчик Р 


для їх зберігання. Важливою особливістю 
задачі є те, що всі контейнери з ВЯП є кон- 
груентними (мають однакову форму та гео- 
метричні розміри) (рис. 3) 1 повинні розта- 
шовуватись лише на одному рівні (не до- 
пускається розміщення контейнерів один 
над одним). 


О А.М. Чугай, С.В. Альохіна 


Рис. 3. Просторова форма 
контейнерів з ВЯП 


Враховуючи вказані особливості, на- 
бір контейнерів К;, іє Їх, можна пред- 
ставити набором конгруентних кругів С; з 
радіусами ї г,іє Їх . Позначимо через 
и; З Со, у;) - координати центру С;. Тоді 


розташування всіх С;, іє Гл, В Е? можна 


визначити за допомогою вектора 


из (цим) єк?М, Надалі круг С., 
трансльований на вектор и; 2 (Х., У;), 
позначимо як С; (и;). 


Технологічні умови розміщення зада- 
ного набору об'єктів в межах зони розмі- 
щення. Такі обмеження можна поділити на 
два типи. 

Перший тип обмежень обумовлений 
забезпеченням необхідних параметрів без- 
пеки, а саме, неперевищення загального 
рівня іонізуючого випромінювання та до- 
тримання режимів теплового зберігання 
відпрацьованого ядерного палива. Для 
опису першого типу умов розміщення 
контейнерів формалізуємо вплив кожного 
контейнеру на загальний тепловий режим 
та рівень іонізуючого випромінення май- 
данчика. Оскільки кожен контейнер ха- 
рактеризується своїми фізичними власти- 
востями (температурою та рівнем іонізую- 
чого випромінювання), які необхідно вра- 
ховувати для дотримання умов безпечного 
зберігання, то кожному С; поставимо у 
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відповідність інтегральний показник К; 
теплової та радіаційної властивостей кон- 
тейнера з ВЯП. Цей параметр визначатиме 
інтегральний рівень впливу кожного кон- 
тейнера на температуру контейнерів, що 
розміщені поруч і, разом з тим, на загаль- 
ний тепловий рівень та рівень іонізуючого 
випромінювання майданичка Р. 
Припускаємо, що значення інтег- 
ральних коефіцієнтів К; визначається екс- 


пертним шляхом. Тоді, відповідно до зна- 
чень К;, іє Їу, розподілимо кола С., 


і є Їм , по группах С;, / Є І,. Нехай кож- 
на група складається з 4;, ієЇ,, кругів, 


де 2 - кількість отриманих груп. 


Для того, щоб мінімізувати взаємний 
вплив іонізуючого випромінювання та 
температурного режиму контейнерів, за- 
дамо мінімально допустимі відстані 


9 . . . -. 
Фі ,,/єЇ г Між колами кожної групи С, 


іє І я» ТА між колами всередині однієї 


групи. 

Таким чином, врахування першого 
типу технологічних обмежень дозволить 
забезпечити підвищення температур па- 
лива та рівня іонізуючого випромінюван- 
ня від конейнерів з ВЯП та забезпечити 
рівномірний розподіл іонізуючого випро- 
мінювання в межах майданича при збері- 
ганні на ньому відпрацьованого палива. 

Другий тип обмежень обумовлений 
забезпеченням умов обслуговування кон- 
тейнерів. Необхідно забезпечити можли- 
вість підходу до кожного контейнера спе- 
ціальної обслуговуючої техніки для того, 
щоб повернути контейнер або перемістити 
контейнер в межах майданчика (рис. 4). 
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поь паса РИчНУ 


Рис. 4. Забезпечення умов обслуговування 
контейнерів транспортером 


Для забезпечення цієї умови необ- 
хідно враховувати розміщення на майдан- 
чику так званої «мережі обслуговування» 
для пересування техніки. Нехай для пере- 
міщення техніки необхідно забезпечити 
наявність смуг з шириною а. Тоді при 
розміщенні об'єктів необхідно забезпечи- 
ти умову торкання С., і Є Гу , смуги зада- 


ної ширини а, що забезпечить підхід 
спеціальної обслуговуючої техніки до 
контейнера. 

Критерій оптимізації розміщення 
об'єктів. У якості критерія оптимізації 
оберемо знаходження максимального за- 
повнення обраного майданчика кругами 
С; з множини І у. 

Таким чином, після формалізації в 
аналітичному вигляді всіх умов оптиміза- 
ційної задачі геометричного проєктуван- 
ня, сформулюємо постановку задачі в на- 
ступному вигляді. 


Задача. Визначити вектор 


2 о 
и з (ці йо угйу) Є В", який гарантує роз- 
міщення максимального числа кругів з на- 
бору С., іє І , у заданій області Р із за- 


безпеченням виконання заданих техноло- 
гічних обмежень (рис. 5.). 
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ай | - зни 
ана нн мм 


Рис. 5. Постановка задачі розміщення 
контейнерів з ВЯП 


Відповідно до постановки задачі, по- 
будуємо математичну модель задачі та ви- 
конаємо дослідження її властивостей. 

Математична модель задачі та її 
властивості 

Однією з найважливіших та най- 
складніших задач комп'ютерного та мате- 
матичного моделювання розглянутого 
класу задач є аналітичний опис взаємодії 
між С; та областю Р. У цій роботі ско- 
ристаємось методом рбі-функцій, пред- 
ставленим, наприклад, в роботах |9, 10). 
Цей метод визнаний на сьогоднішній день 
найефективнішим для розв'язання подіб- 
них задач. 

З метою формалізації умов розмі- 
щення С; в межах Р за допомогою ррі- 
функцій побудуємо наступну множину 


О А.М. Чугай, С.В. Альохіна 


Са св? В). Ця множина завжди може 


бути представлена кінцевим об'єднанням 
первинних об'єктів о , і 21,2,3,4, тобто 


б у 5 
са.|)0, 1 1)95, 1 1)9з; 
дз! чі З 
С 
|) ом 
а 


де 

Оз - (о у)є В? г ду (5 у) 20), 

і є 15 « (1,2,...,6), 

діб У) загух ЯБузУ КС), 

Оз; з (о у) є В? 1 да я (5 у) 20,1 21,2), 
іє 7, - 11,2,..5 У 

Хо 5 У) з Фр Нв ЛУ ЧОСо й, 

Озу «(5 у) є В бах, у) 20), 
зерно» 

аву) втї; «(кох ЗОЗ У, 
Олі з (Со у)є В г дда (ху) 20, 

од (х, у) 20,1 - 1,2), 

ер цій 

Ха 05 У) зад На ДУ НОД 
сдібьу (кох ОЗ) їі» 
Р(А,Ву22к, р(А,В) - це відстань між 
точками А і В, тобто 0, 6 напівплощи- 


ною, 0»; - опуклий конус; Оз; - коло і 


Од і - перетин напівплощини та допов- 


нення кола до К/ (рис. 6). 
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Хал(Х, у)-0 


ФС 7-0 ря? у)-0 


Рис. 6. Об'єкти 0,,,0,5 0. О,, 


Наприклад, набір С для області Р,, 


зображений на рис. 1, може бути отрима- 
ний у такий спосіб (рис. 7): 


е-|оа бе Мо Ме) 


Рис. 7. Представлення множини С 


Для розв'язання задачі пропонуємо 
підхід, який дозволяє звести розв'язання 
задачі до розв'язання послідовності задач 
з лінійними цільовими функціями. Для 
цього радіуси г з С;, іє І, приймаються 
змінними. Таким чином, радіуси форму- 
ють вектор у" «(п,р,азг)Єк". 

Значеннями рбі-функції є деякі міри 
як перетину двох геометричних об'єктів, 
так 1 найкоротшої відстані між ними, за- 
лежно від їх взаємного розташування в 
просторі. 
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Використовуючи  ррі-функції для 
первинних об'єктів та складних 2Д-об'єк- 
тів, математичну модель послідовності за- 
дач можна представити наступним чином: 


п 
ЕХ" уз тах КОХ") с тах У п, (1) 
і-І 
Х" «(и",у")є М, 
п -1,2,....741, 


і) о ЗУлА ска 0) 
Фу(йрй з Гі) 2 0,1, є Іс, 
г-г20,в 20,ієЇу, 
де 
2 2 
Фу(цьньтьТі) 2 (0 хо) Оу -Ур- 
8 92 й 
во іні? 
Фіи:г)  -  рбі-функція для С; і 
С, «сІ(ВОХРУ, 
В гло ол ло - ос 2 - с а2 
Фу(й Й з, їі, Гі) -( -Хі) РК -У,) - 
«(бабу а), 
Іс - множина кластерів, які включають 
кола С;. Функція Ф(и;, г) описує прина- 
лежність С; до області Р. 
Ррі-функцію для С; 1 С -сі(К? Р) 
можна представити наступним чином: 
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Ф(иї, п)з тіп (ФІ; і ЗСя г), 
ФС ціп), 1є1 Я 
де 


Мп) піп ФОС (шіьв), 


и (и, п), зеіучєт| 
зі Чи г) - ррі-функція для С; 1 А; 
Фі (шоп) (ох Ого У) з 
- зр 
Фіх (ці, г) - ріі-функція для С; і СІ,; 


- М (цьв)- отах 
іа ай кіт 


(тах (тіп У ці (и, г) 


Фидк (і Ов, РОЛ Ота о 
Фі зн СР г) є рбі-функцією для С; 1 Мі,; 


ФО 


(ші) - піт |ФУО Чи, г), ) є Їз, 


зн зва Фе 0, 
3 а и,г), і 
дО г) є ррі-функцією дляС; 1С; 
порно дат) є-д цу) 


пит в) тах (піп (/у (шт г), Фу(иі г) 


до, п); Хо із ці трі 


ФіЗ(иьпуа (ох Оу з 
«Рі , 
ф- (цьг)з тах | фу (щу ту), фуз(ці» ті), 


ння ТЕТОЛА г), 


жалі п), ду иьті) 


О А.М. Чугай, С.В. Альохіна 


Фі" цуп), ФЕЗЗ(цд,п), 
Фі цу), Фі ит г) мо ріші- функції 
для С; і 0; С б 02; С в ОС Чі 
Од) відповідно. 


Компоненти функції Фі; и, г) 


мають вигляд: 

фФіїиь г) - тіп (с, ал (иг 1), Фігз (й, Г), 
фол (пу (пз) (ші) 

Фіз (шьту) є поїт| Фут ці» пі) фозд(ші» ті); 
фузщьті) є тіп (аузз (шо ті) Муз г); 
олшьпує (С охаї ЧОрЗУЛУ ЗТ 


2 2 
Фіз(ирт) 2 (ох) НКУ; З У) -п 


аззцьг) З (бу «бу ЗО ох Зо) 
ул) задх; лу; КСД, 
ау З Уа 7 Уїз» Раз "Суд Ху 
а зх НВ з У;3)» 
уна(нь г) зах нвоур кс, 
аз» 7 Уу6 7 Уув» фуд З "(Ху6 7 Ху9), 
сі - азу -Ві2Удз5), 
Міз(ці п) га зх; КБУ; КС, 
азз 7 Уу5 7 Уа» Фіз З "(Х)5 3-Х), 
сіз З "(аззхо -Фузує), 
ууд(цьг) зад КВУ; НСД, 
а 7 Уу3  Уїт» Рід 3 "(23 7 Х)т), 
с4зЧа хз ВУ) 
РЕТОТИ г)уза із нФуУр КО у 
аі5 З У;8 7 У ув» Ф)5 З З(Х)8 7 Х16)» 


сів З -(аз9Хі6 "РузУ 6)» 
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даліц бує -Далц)-т, 
дадіиьн) є "Дазаці)-п, 
(хльУп (хід Уго)» (хз, Ура)» 
(х ід» Уза)» (Х 15, Уі5)» (хув» Уу6)» 
(хіть У; 1 (Х)в» Уув) - координати точок 
А, В, В, Е,Е, С, М та М, відповідно. 
Тре 0-Іеуеі ої ФІЗ (и,) 


Тре 0-Іеуеі ої Фе (и,) 


Рис. 8. 0-рівень Фе (и,) 
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Якщо границя множини Од; утво- 
рюється однією прямою лінією, то 


с ; 
ФО (и;) можна написати простіше 


С 
ФІ (ціп) тах | фу (мі г), 
Й 5 
Фіз(иї; ті) ал, '?) 
Ф и: г)20, якщо хоча б одна із систем 
нерівності виду 
С . 
Фе ирп)20,/є 15, 


С ; 
Фе иьпу2 іє, 


со | 
Ф"(и,,г)20, є хх, 
ГУ гра і) ї»П / Є 


1? 


Фет) 20, / є /;, 


1211 
СС 
Фіїг (иьп)201ЇєЇ,55 єї» 
СМ 
Фіїд (ц,п)201є1,,4є Ла, 
ФО 20 є аб 
виконується, де Фі" (и,п)20 є або 
. й хі ж 
однією з нерівностей | До д(и, )20, 
ж . - 
Жіз у п)20 або система нерівностей 
ау (ць) 20, 


Фа пу20 є  ц або 


в однією з 
нерівностей  Фуз(и;,п)20, 
далі п)20, дуг) 20 або одна 
(и) 20, 
Фіз (ні п) 20, 
із систем нерівностей 4 (и;,)20, 
Міг (ці, 5) 20, 
Мііз(щьт) 20, 
ОСТ РА Базрнн 
ни уз (іт) 20; 
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Фі з ,0)20 або нерівність, 
Ра ідк (5) 20 або система нерівностей 
Фідк (пі) 20, 


Мах (б п)20. 


Властивості математичної моделі 

Виконаємо дослідження деяких важ- 
ливих властивостей математичної моделі 
задачі 


п 
М 
ЕХ" )етах Е(Х") зтах У г , 0) 
і-1 
Х" «и",у")є М, 
Якщо Е(Х" )-пг, тоді точка Х" 
є глобальним максимумом задачі. 
1. Якщо м СХОВУ «(бук і 
Ж ж 
ХУ зу еру де Ж ОСУ" 
і Х"" є глобальними максимумами за- 
дачі (3), тоді розв'язання поставленої 


. . ож ж 
задачі досягається в точці и" зи" , яко- 
му відповідає г розміщених кругів. 


2. Оскільки Е,(Х") є лінійною, то ло- 
кальні максимуми досягаються в край- 
ніх точках У/,,. 


3. У/, задається лінійними та нелінійни- 


ми нерівностями. 
4. Ф(и;,п)20, якщо хоча б одна із сис- 


тем нерівності виду 


Фі іп п)20, є я, 
ФІипу20,/є у, 


ФЕЗ(шп) 20, 
- є є, 
Гшупуз б і 


С 
про СРО 
й Си, г)201єЇ,,8 єї, 


Фій М г)20,1єІ,,4є 13, 


СО» (п; 


виконується, де Фі; пу20 є або 


їз Й 


однією з нерівностей жо л(цьг)у20, 


О А.М. Чугай, С.В. Альохіна 


; 5 г; 
Жіо іди п)20 або система нерівностей 


оіщуп)20, збо 
мі (5 п)20; 


однією з нерівностей 


гву20 є або 


ї» й 


Фаз(цьп)20, 
ж ж 

дали пу20, діди, п)20 або одна 
ол (ці г) 20, 


із (ціз їі ) 2.0, 


із систем нерівностей 4 (и; 1)20, 
М ії2 (ці; , г) 20, 
уз (щі, п)г 0, 


оіт(ць г) 20, Фіз (и, п) 20, 
ни п)20, уз (ць п) 20; 


М и; ,г)аг0 або нерівність 


Фій 

Ж (и: п)20 або система нерівностей 
Фк (Щі п) 20, 
Мак (Мі г) 20. 


1. З попереднього пункту випливає, що 
Фіи:,п) включає функцію 


ФІ ; 
пої» леді, 
ФІ 
й я 22)» / -і, ді з 
Фі зав ОНИ ву, Феї 2 91283 
кожна з яких є функцією виду 
С 
Фі Чи пу) я тах (фу (цу пу)фул (шт), 
Фіз(иі, 1), Ф ТАГОТИ їі), 
ж к 
діал діауа гі) 
Це означає, що якщо хоча б одна з 


функцій у фігурних дужках є позитивною, 
тоді відповідна функція 


Со Со 
Фі (ціп). Фі ""(ці,п) або ФІ (ці, п) 
є позитивною. Тому Ф(и;п)20, якщо 
виконується хоча б одна із систем 
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нерівностей Гі(ип)20, 


с 9 
Ея ні: 
5 є Ід - (12, 237-69-ЩЦИту). 
ї-1 4-1 
2. На основі попереднього пункту та ма- 
тематичної моделі (3) будь-яке розмі- 
щення п кругів у Р описується систе- 
мою нерівностей виду 
ГА (ць) 20 ре. 
Фу(цьнітьті)20, і,7єГ,, іх), 


перо реїоз 


в2б, 1єЇ,. 
Система визначає певну підобласть 
У | Кожна система нерівностей 
складається не менше, ніж З 


п(Ачси чн) ЧО -3)) лінійних / та 


нелінійних нерівностей. 
3. За пунктом 5, кількість підобластей 


становить 7-7. Таким чином, 


область У/, можна представити як 


7 
У, з |) ШИЮ 
т-І 


Слід зазначити, що серед 7) систем 
існує ряд несумісних систем. 


4. Якщо точка Х"! є Міо г Локальний 


максимум щодо У/,; і одночасно 


ід 
хі" ЄМ; , /є Л с Ц1,2,...По| то не- 


з ж М 
обхідно довести, що Х'"" локальний 
максимум щодо У), . 


5. Задача (1)-(2) є КР важкою. 

Висновки 

У роботі задача оптимального розмі- 
щення контейнерів з ВЯП з урахуванням 
заданих технологічних обмежень сформу- 
льована у вигляді задачі оптимізації гео- 
метричного проєктування. Виконано опис 
всіх умов задачі геометричного 
проєктування. 

Побудовано математичну модель за- 
дачі упаковки конгруентних кругів у бага- 
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тозв'язну область, межа якої складається з 
дуг кіл та відрізків ліній. Математична мо- 
дель представлена у вигляді нелінійної за- 
дачі оптимізації. 

Метод рбі-функцій дозволяє нам по- 
будувати математичну модель проблеми, 
в якій область припустимих розв'язків мо- 
же бути представлена у вигляді об'єднан- 
ня підобластей. Кожна з підобластей опи- 
сується системами нерівностей, ліві части- 
ни яких є неперервними функціями. Це 
дозволяє використовувати сучасні методи 
нелінійного програмування для розв'язан- 
ня задачі. 
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